On the heat current in the magnetic field: 
Nernst-Ettingshausen effect above the superconducting transition 
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For maintaining gauge invariance in a magnetic field, the heat current operator should include the 
magnetic term. Taking this term into account, we revised calculations of the Nernst-Ettingshausen 
effect above the superconducting transition. We found that the fluctuations of the modulus of the 
order parameter do not change the particle-hole asymmetry (PHA) of the thermomagnetic effects. 
As in the normal state, the thermomagnetic effects in the fluctuation region are proportional to the 
square of PHA and, therefore, small. Large Nernst effect observed in the high- Tc cuprates requires 
vortex-like excitations due to the phase fluctuations, which are beyond the Gaussian-fluctuation 
theory. 



PACS numbers: PACS numbers: 71.10.-w 

The observation of large thermomagnetic effects in 
cuprates turns out to be crucial for understanding high- 
Tc superconductivity In the Nernst effect, the 

electric current is induced by crossed temperature gra- 
dient and magnetic field, 'f = N\SIT x H]. In an open 
circuit, the Nernst voltage Ey is given by 



N. 



(1) 



xy 



where r],j;x is the thermoelectric coefficient, and a. 
are the electrical and Hall conductivities. 

In metals, the Nernst voltage is small, because it re- 
quires strong particle-hole asymmetry (PHA) in the elec- 
tron systems. PHA of the thermoelectric coefficient rj^x 
is given by the factor T/ep, where is the Fermi en- 
ergy; PHA of the Hall conductivity axy may be presented 
as r2//r, where ^Ih ~ eH/mc and r is the electron mo- 
mentum relaxation time. Coefficient Nxy contains both 
PHA factors of r]xx and axy Thus, in Eq. ^both terms 
are proportional to the square of PHA. Moreover, if r is 
independent on energy, they cancel each other. 

The important result of Ref. is that the large mea- 
sured Nernst voltage originates from Nxy, i.e. the first 
term in Eq. ^ strongly dominates over the second term 
in the temperature range extending from below Tc to a 
broad interval above. It now seems clear that the re- 
gion above Tc is in fact the vortex fluid phase. Various 
phenomenological models are suggested to describe this 
phase and heat transfer by vortices 0, IE IE • It is also 
believed (see a review 0), that, at least close to Tc, the 
Nernst voltage is described by the microscopic Gaussian- 
fluctuation theory, which pre dicts the finite Nxy in the 
zero order in PHA |E IE ^E Below we will show 

that this statement is wrong. In the Gaussian model 



with ordinary particle-hole excitations, it is not possible 
to change the PHA of Nxy. To obtain the large Nernst 
effect due to vortex-like excitations, it is necessary to go 
beyond fluctuations of the modulus of the order param- 
eter. 

The Gaussian-fluctuation theory allows one to calcu- 
late the Ettingshausen coefficient T, which describes the 
heat flux induced by crossed electric and magnetic fields, 
j/i _ — X[E X H]. According to the Onsager relation, 

N = T/T (2) 

Results of the previous theoretical works can be for- 
mulated by two statements. First, in the infinite sample 
the coefficient T near Tc does not require PHA, 



ST. 



inf 
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for 3D, 



(3) 



where r? = (T - T,)/T, and a = {T - Tc)UT)Y, ^{T) is 
the coherence length. For the first time, Eq. |2|has been 
obtained in Ref. |S| using the time-dependent Ginzburg- 
Landau (TDGL) theory. Later it has been reproduced in 
Ref. ^E- Finally, Eq. |3|has been confirmed by the di- 
agrammatic calculations in the frame of the Aslamazov- 
Larkin (AL) approximation 

Second, according to Refs. 0, Q, 0| in the finite 
sample the Eq. 13 should be corrected due to surface 
currents related to thermal magnetization currents, 



ST = 




for 2D, 
for 3D. 



(4) 
(5) 



where the /i is the fluctuation magnetic susceptibility 
[l3j |. This correction would decrease the T to one third 
of its value given by Eq. El 
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In the current paper we develop consistent microscopic 
description of the heat current in the magnetic field. Cal- 
culating by the Kubo method and SN by the trans- 
port equation, we show that both statements expressed 
by Eqs. |31 and 01 are wrong. Our conclusion is that in 
the frame of Gaussian-fluctuation theory the coefficients 
ST = SN = in the the zero order in PHA. 

We start our calculations with discussion of the energy 
and heat current operators. The energy current of in- 
teracting electrons may be obtained from the equation of 
motion for the electron field operators or from the energy- 
momentum tensor 0|. For electrons interacting in the 
Cooper channel, the energy and heat current operators 
have been derived in our work In the external fields 
H = i[]<.H X Ah] and E = — ik^;^, the electron energy 
has a gauge-invariant form 



= (p + eAH/cf/2m + e(f>. 



(6) 



The energy current operator is easily generalized to in- 
clude the external fields. 



J' = X! '^^P^p"? + — (vAH)a+ap + ve(j)a+ap 
p p p 

+ ^ v/x a+Qp - A/2 ^ (v -t- v') a++p,_p»a+,ap.ap 
p p,p':P" 
+ (v + v')f/.™p(R») exp[i(p-t-p')i?»] a+flp' , (7) 

p,p',Ri 

where and Cp are the electron creation and annihila- 
tion operators, — p'^/2m ~ ^ is the chemical poten- 
tial), Uimp is the impurity potential. The first and the 
forth terms in the Eq. ddescribes the energy flux of non- 
interacting electrons without fields. The second and third 
terms are due to electron interaction with magnetic and 
electric fields. Two last terms are many-body corrections 
due to the electron-electron and electron-impurity inter- 
ations. According to Ref. 0, these two terms generate 
diagram blocks (AL-blocks) proportional to PHA, but 
contributions of all blocks with PHA cancel each other. 
Anyway, these many-body corrections can be neglected, 
if T is calculated in the zero order in PHA. 

The thermal energy can be defined as the electron en- 
ergy counted from the electro-chemical potential e(f) + fi. 
Then the heat current may be presented in terms of the 
energy and electric currents as J'' = J'^ — {fi/e + 0)Je. 
Obviously, the third and the forth terms in Eq. do 
not contribute to the heat current. Thus, calculating the 
heat current we should take into account only two heat 
current vertices 7^' and 72 corresponding to the first (ki- 
netic) and second (magnetic) terms in Eq. [3 

In the linear response method, the thermoelectric co- 
efficient is given by the correlator of the heat and charge 
currents Calculating response to E x H, it is conve- 
nient to use the gauge conditions -Ah — k^; • A// = 0, 
in this case E x H = —iAhO^h ■ E). 



Two leading terms in the heat current operator gen- 
erate two diagrams, which describe T in the Aslamazov- 
Larkin (AL) approximation. The left blocks and Bj 
in Figs, l.a and l.b are blocks of electron Green func- 
tions connected with heat current vertices 7^* (kinetic) 
and 72 (magnetic). General for both diagrams, the right 
block B*^ includes the electric current vertex 7"^ — ev-e^;, 
e^; — 'Ei/E. Block B^ includes the magnetic vertex 
7^ = (e/c)v • Ah- The wavy lines correspond to the 
fluctuation propagator, which is given by 



respond 



L 



R.A 



(g,a;) = (A-i-F«'^(g,c^)) \ 



(8) 



where A is a constant of the interaction in the Cooper 
channel and P{q, lo) is the polarization operator. 



2 \ ttT 



- a?" ± ^ + 7^ ) , (9) 



V is the electron density of states, is the Debye fre- 
quency, and is the Euler constant. The last term in 
Eq. Ois proportional to PHA and can be ignored in zero 
order in PHA calculations. 

For an arbitrary electron momentum relaxation time t, 
the AL blocks B'^'''-^ of electron Green functions G^'^^^ 
with vertices 7^^, 7^, and are given by [l^ll5| 



-rie,h,H T I dp de 
rJ,- = im / — — — V 7, 



(2^)^ 



e,H,H (10) 



(1-0^ 



5(e) = -tanh— , C - — 



dp 



TTi/r J (2^)3 ^^P^i'P' ^ > 
[G^r ^{e-S,, + il2T)-\ (12) 



The block B'^ with the electric current vertex, 7*^ = ev • 
e^, may be presented as 0,0] 

B^(q) = 2e VqP^(q,0) ■ GB = 2ei/aq-eE. (13) 

The block B-^ with the vertex 7-^ — (e/c)v • Ah is 

B^(q) = 2{e/c)vaci-AH. (14) 

The block B^* with the kinetic heat current vertex, 



• e^f, (Gj-h = li^ ||A), is given by jl5j (see also [lOj) 

'S,\{<\,ijj) = uj Vt^P'^{q,0) ■ e^h = ujiya q- e^h. (15) 

Next, we calculate the block Bj with the magnetic heat 
current vertex 72 = (v • A/f)(v • Sjh). The integral over 
angles of the electron momentum involves only the vertex 
72 , because the heat current is in the direction of Ah- 
To obtain an imaginary part in Eq. 1101 the integral 



d^ {Gp )^Gq_p 



2iTi 



(2e — — q • V — i/tY ' 



(16) 



should be expanded in lo (in calculations of Bg it is ex- 
panded in q • v) . Finally, we get 

B^'(q,c^) = 2Ah - VqP^(q,0) = 2{e/c)iuvaAH- (17) 
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The first AL diagram (Fig. l.a) for the Ettingshausen 
coefficient is based on the blocks and B*^, the vector 
potential A enters the fluctuation propagator. This di- 
agram was calculated in Ref. |lll | and its contribution 
to Tinr is given by Eq. O The same result has been 
obtained in the TDGL formalism in Refs. [8lllo|. 

The second AL diagram (Fig. l.b) is based on the 
blocks B2 and B*^. The analytical expression correspond- 
ing to this diagram is 



y(2) TJ _ Ck 



dd du: ^2^\L<iLi^LlL% (18) 



(27r)" 2ti 2n 



where = coth(w/2r)(L^ - i^), L± is used for 
L(q ± k/2, u! ± ft/2), and n is the system dimensionality 
with respect to the coherence length ^(T). Expanding 
the integrant to the linear order in f2 and k and calcu- 
lating the integrals over uj and q, we find that the con- 
tribution of the second diagram, T -^^^ cancels completely 
the contribution of the first one (Eq. ^ . The analogous 
cancellation of two diagrams with the heat current ver- 
tices and 72 takes place in the case of noninteracting 
electrons (see Appendix in Ref. [l^). 

Thus, without PHA the Ettingshausen effect is absent, 
T = 0. To get nonzero result, wc should expand the 
fluctuation propagator (Eq. (SJ up to the second order in 
PHA (in the ffi'st order in PHA both diagrams give zero 
results), i.e. we expand the polarization operator (Eq. O 
to the second order in jlu. In the second order in PHA, 
the thermomagnetic coefficients are 



T 47rc V TT y \ a/o/^ for 3D; ^ ' 



1 aini 



In ^"^^^ [13. Thus, the Etting- 



where 7 = 

2eF omeF ttTc 
shausen and Nernst coefficients in the fluctuation region 

are proportional to (T/cf)^- Taking into account that 
in this region the thermoelectric coefficient ij^x and the 
Hall conductivity axy are proportional to {T/ep), we see 
that in Eq. for the Nernst voltage both terms are of 
the same order, {T/ep)^. 

The above results obtained for the infinite sample (or 
in the interior of the sample) do not show any con- 
tribution of the magnetization currents and, therefore, 
in the finite sample no additional corrections are re- 
quired. Because the role of magnetization currents is 
widely misinterpreted in literature, below we briefly an- 
alyze the energy current induced by the magnetization 
current and show that this term is erroneously attributed 
to the heat flux (Eq. In the potential (/'(r), the elec- 
tric magnetization curren t transfers the energy flux 



imag 
'f 

•imag 



4>3mag (Eq- 37 in [I8j). Taking into account that 



ciik Ah, we get 



Jmag 



c^[H X E]. 



(20) 



In the linear response method, this term is given by the 
diagram with the energy current vertex 7*^ = e</) v (the 
third term in the energy current operator, Eq. [TJ and 
the magnetic vertex 7^ = {e/c){AH ■ v). Without (j) this 
diagram gives the jf^^g, and, therefore, we get j^^g = 

0j 



mag ' 



As we discussed, the vertex 7* = evcj) does not 
contribute to the heat current, because the heat energy 
is counted from the electro-chemical potential. 

While the magnetization currents do not contribute 
to the heat transfer, they play an important role in the 
charge transfer. In the interior of the sample the elec- 
tric current consists of the transport and magnetization 
components, j^„^ = j^^ + imag- According to Ref. |l8| . 



Jmag 



(VT X H). 



(21) 



The magnetization currents are divergence-free. There- 
fore, the total magnetization current through the sample 
cross-section must be zero, i.e. the bulk magnetization 
currents are canceled by the surface currents. Therefore, 
in the finite sample, the measured Nernst coefficient N is 
determined by the transport currents, N = jf^/[Vr x H] 
. Calculated from the theory, the Nernst coefficient in 
the infinite sample may be associated with the current in 
the interior of the finite sample, iVm/ = Jinr/I^^ ^ H]. 
Using Eq. [21 we get 



SN = 



Jmag 



[VT X H] 



ONinf - C — 



dT' 



(22) 



The coefficient SNmf may be calculated by the quan- 
tum transport equation, where VT can be easy incor- 
porated [ij, |l^ |l^ llSl |2(]j (as a nonmechanical pertur- 
bation, VT cannot be incorporated into the Kubo for- 
malism). In the AL approximation, calculation of the 
Nernst coefficient is analogous to calculation of the Hall 
effect with the only difference that the derivatives 
of the fluctuation propagator L^'^{q,uj) with respect to 
Lu in the Hall effect should be replaced by the temper- 
ature derivative in calculations of SNi^f- Taking into 
account temperature dependencies of all coefficients in 
the fluctuation propagator (Egs. |Hland|5I, for the 2D- 
superconductor we get 



6N,, 



inf 



SN- 



e I a 
3ttc \ rf' 



a 



a da 



(23) 



Thus, in the inflnite sample or in the interior of the finite 
sample, the Nernst coefficient consists of a large term 
(the second term in Eq. I23|) . which is nonzero to zero 
order in PHA. However, in the finite sample this term is 
canceled by the contribution of the magnetization cur- 
rents, cdn/dT. According to Eq. the rest is equal to 
SN, which satisfies to the Onsager relation SN = ST/T. 

Thus, in the finite sample both Nerst and Etting- 
shausen coefficients in the fluctuation region (Eq. I19|) 
as well as the coefficients in the normal state {Nq and 
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To) are proportional to the square of PHA. Therefore, 
the relative corrections do not consist of any PHA fac- 
tors. For example, for 2D superconductors we obtain 

^ ^ ^ ^ i I (24) 

TVo To -n niax{T,l/r}- ^ ^ 

These corrections are of the same order as the relative 
corrections to the conductivity, themoelectric and Hall 
coefficients. Our consideration was limited by the one- 
band model. However, Eq. |31 is also valid for the 
two-band model, where the electron and hole contribu- 
tions are additive and thermomagnetic effects are large. 
In other words, our key statement is that the interelec- 
tron interaction cannot change PHA requirements for the 
thermomagnetic coefficients. The previous works gave 
5N/Nq with an additional large factor of (e^/T)^. How- 
ever, none of experiments has shown such huge effect in 
ordinary superconductors. While claiming opposite, the 
very recent work on SiN actually supports our re- 
sults in the form of Eq. |31 if the authors would have 
normalized SN obtained in the fluctuation region by A^o 
taken from measurements far from the transition, say 
at T ^ 2Tc. In fact, they evaluated A'o by the term 
Vxx'^xy/'^xx, which was found to be many orders of mag- 
nitude smaller than the measured N even far above Tc- 

Theoretical problems of thermomagnetic transport 
rests on identification of the correct form of the heat cur- 
rent operator, which takes into account the interaction 
corrections and corrections due to external fields. As 
we have seen above, the transfer of the potential energy 
e(/)(r) by the magnetization currents contribute to the 
energy current, but not to the heat current, because the 
electron thermal energy is counted from fi+e4>. Thus, the 
magnetization currents do not affect the Ettingshausen 
effect. On the other hand, the magnetic field gives an im- 
portant contribution to the heat current operator. Tak- 
ing into account the magnetic term - viz., (e/c)(v- A/f)v 
(Eq. [T]), - we have shown that the corresponding correc- 
tion to the Ettingshausen coefficient cancels completely 
the result (Eq. |3l obtained with the kinetic term in the 
heat current operator. Thus, previous calculations erro- 
neously attributed the electric term in the energy current 
to the heat current 0, ^ and overlooked the mag- 
netic term [lISESOl- 

In summary, we have shown that any calculations of 
many-body effects in thermomagnetic phenomena should 
include the magnetic term in the heat current operator 
(the second term in Eq. [7J . Magnetization currents in the 
electric field contribute to the charge and energy transfer, 
but not to the heat current. Only in this way, one can 



obtained the Nernst and Ettingshausen coefficients that 
satisfy to the Onsager relation (see Eqs. 1221 and l^^ . The 
Gaussian fluctuations, i.e. the fluctuations of the mod- 
ulus of the order parameter, do not enhance the ther- 
momagnetic effects, because these fluctuations do not 
change character of elementary excitations (particle-hole) 
and, therefore, do not change the particle-hole asymme- 
try of the thermomagnetic coefficients. The large Nernst 
effect observed in high- Tc cupprates requires vortex-like 
excitations related to the phase fluctuations. 
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Fig. 1. 



Fluctuation AL diagrams describing the heat current in 

crossed electric and magnetic fields. Wavy lines stand for the 
fluctuation propagators and straight lines stand for the electron 
Green functions, which form the AL blocks. 



